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Recently, a simple non-interacting-electron model, combining local quantum tunneling via quan- 
tum point contacts and global classical percolation, has been introduced in order to describe the 
observed "metal-insulator transition" in two dimensions jlj. Here, based upon that model, a two- 
species-percolation scaling theory is introduced and compared to the experimental data. The two 
species in this model are, on one hand, the "metallic" point contacts, whose critical energy lies 
below the Fermi energy, and on the other hand, the insulating quantum point contacts. It is shown 
that many features of the experiments, such as the exponential dependence of the resistance on 
temperature on the metallic side, the linear dependence of the exponent on density, the e 2 /h scale 
of the critical resistance, the quenching of the metallic phase by a parallel magnetic field and the 
non-monotonic dependence of the critical density on a perpendicular magnetic field, can be naturally 
explained by the model. Moreover, details such as the nonmonotonic dependence of the resistance 
on temperature or the inflection point of the resistance vs. parallel magnetic are also a natural 
consequence of the theory. The calculated parallel field dependence of the critical density agrees 
excellently with experiments, and is used to deduce an experimental value of the confining energy in 
the vertical direction. It is also shown that the resistance on the "metallic" side can decrease with 
decreasing temperature by an arbitrary factor in the degenerate regime (T < Ef). 

PACS numbers: 71.30.+h, 73.40. Qv, 73. 50. Jt 
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I. BACKGROUND AND INTRODUCTION OF 
THE MODEL 

The surprising experimental observation of a metal- 
insulator transition in two dimensions [§fj6|, in contra- 
diction with the predictions of single-parameter scaling 
theory for noninteracting electrons JtJ, has been a sub- 
ject of extensive investigation in recent years. Theories 
ranging from attributing the effect to scattering by im- 
purities |sj to those suggesting "a new form of matter" 
H have been proposed. Some theories, based on the 
treatment of disorder and electron-electron interactions 
by Finkelstein [|l0|, have been put forward [jnj, while 
other approaches considered spin-orbit scattering JlJ] or 
percolation of electron-hole liquid jL3|. Altshuler et al. 
H gave several arguments why this transition is not due 
to a non-Fermi liquid behavior, including the fact that 
the exponential increase of the conductance with tem- 
perature persists to high densities where the conductance 
is almost two orders of magnitude larger than the criti- 
cal conductance, and the fact that the Hall resistance is 
rather insensitive to temperature, and does not display 
any critical behavior. Some experimental results sup- 
porting the conclusion that the transition is not driven 
by interactions that were mentioned in |jj included the 
fact that such a transition was observed also in high- 
density electron gas upon the introduction of artificial 
disorder [|l4| and the fact that increasing the density in a 
parallel electron gas increases the conductance ||, even 
though the interactions are screened by the parallel gas. 
More recently, the compressibility on the metallic side 
of the transition was measured (T^j and was shown to 
be accurately described by Hartree-Fock approximation, 



again indicating a normal Fermi-liquid behavior. Sev- 
eral other recent experiments ]l6| , pTt have demonstrated 
weak-localization behavior on the metallic side with very 
little effect of electron-electron interactions. 



Recently Q I proposed a simple non-interacting elec- 
tron model, combining local quantum tunneling and 
global classical percolation, to explain several features 
of the experimental observations. At low electron or 
hole densities the potential fluctuations due to the disor- 
der cannot be screened and they define density puddles 
(density separation into puddles in gated GaAs was in- 
deed observed experimentally by Eytan et al. p8|, using 
near-field spectroscopy). These puddles are connected 
via saddle points, or quantum points contacts (QPCs). 
It is now established that even at low temperatures and 
for open puddles (or quantum dots), the dephasing time 
may be shorter than the escape time from the puddle 
p9| . Thus it is assumed that between tunneling events 
through the QPCs dephasing takes place, and the con- 
ductance of the system will be determined by adding clas- 
sically these quantum resistors. (A related model was 
introduced by Shimshoni et al. to describe success- 
fully transport in the quantum Hall regime.) Each saddle 
point is characterized by its critical energy e c , such that 
the transmission through it is given by T(e) = 0(e — e c ). 
(I assume that the energy scale over which the transmis- 
sion changes from zero to unity is smaller than the other 
relevant energy scales, to avoid additional parameters). 
Then the conductance through a QPC is given by the 
Landaucr formula, 
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where fj, is the chemical potential, and /fd is the Fermi- 
Dirac distribution function. 

The system is now composed of classical resistors, 
where the resistance of each one of them is given by (Q) , 
with random QPC energies. 



so that in the case Ui — > (a regular random resis- 
tor network), a(n) ~ & m ( n ~ n cf ; , while in the case 
<7 m — > oo (a mixture of an insulator and a supercon- 
ductor), a(n) ~ (Ji(n c — n)~ t . (In the above I used the 
notation x l = sign^)^ 1 .) The exact form of f(x) is 
not very important, and in the following I have chosen 
f(x) = log(B + exp[a;])/log(B + exp[-a:]), with B = 2. 



III. "ZERO" TEMPERATURE 



II. TWO-SPECIES SCALING THEORY 

In [Q I presented numerical calculations to be com- 
pared to the experimental data. Here I present a dif- 
ferent approach, based on the scaling theory of a two- 
species percolation network. At low temperatures the 
resistors can be divided into two groups, the conduct- 
ing ones (e c < /i), whose conductance is about 2e 2 /h, 
and the insulating ones (e c > fi), whose conductance is 
nearly zero. Thus the distribution of the conductances 
will be a two-peak distribution, where the weight of each 
peak will be determined mainly by density (or chemi- 
cal potential) and the position of the conducting peaks 
will be determined mainly by temperature. Since most 
properties of such a percolating network are insensitive 
to details of this distribution, I replace it by a two-delta 
function distribution, namely replace the network by a 
network comprising of two types of conductors: 
An effective conductor a m describing the metallic QPCs, 
and whose conductance is given by (Q), with an appro- 
priately averaged e c : 

2e 2 1 
Gm ~ h l + exp[-A/fcT] [) 

(A, which depends on the potential fluctuations distri- 
bution, is taken as unity in the following, i.e. it defines 
the temperature scale) , and an effective conductor <7j de- 
scribing the contribution of the insulating phase. The 
conductance of the insulating QPC is dominated by acti- 
vation PHI], cj{ — a a exp[— A\/T]. (Indeed, experimental 
invesigations reported that "two different contributions 
to the conductivity (or two conducting systems) may ex- 
ist, one with a metallic temperature behavior and an- 
other one with a standard, insulating, weak-localization 
behavior" (l7).) 

The scaling form of the two-dimensional conductance 
of such a two-phase mixture, a, near the percolation 
threshold is well known I221J2T 



(3) 



with t ~ 1.3, the conductance critical exponent for two- 
dimensional percolation, and 



At low enough temperatures, such that T <C A, A\ 
the conductors have either zero conductance or a con- 
ductance equal to 2e 2 /h. If the dephasing time is still 
finite at these temperatures, one has a random-resistor 
network, which exhibits a second-order percolation tran- 
sition p3| |. In Fig. 1 I fit the lowest- temperature experi- 
mental data of (|-[J to the expected critical dependence. 
Clearly, the agreement with the classical percolation pre- 
diction is excellent. Such as agreement with the classical 
percolation critical behavior may serve as an additional 
experimental indication that the dephasing is still finite 
at the lowest available experimental temperatures. 
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Fig. 1. Comparison of the lowest temperature data of [g| (two 
sets of data, triangles and squares, 330mK, density given by 
the lower axis) and of ji| (circles, 57mK, density given by the 
upper axis), and of the n-type data || (diamonds) to the pre- 
diction of percolation theory (solid line). Inset: Logarithmic 
derivative of the data [[| which gives a line whose slope is 
inverse of the critical exponent. The percolation prediction 
(t ~ 1.3) is given by the solid line. For comparison a t — 1 
slope is also shown (broken line). 
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In the inset I plot the experimental data 



for 



l/(<ilog(7/<ilogn), which, if indeed a ~ (n— n c )', is given 
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by (n— n c ) /t. The data indeed fits on a straight line, with 
a slope given by 1/t ~ 1/1.3. For comparison a straight 
line with a slope of unity is also depicted, in order to 
demonstrate that a critical exponent of unity cannot fit 
the data. 



IV. TEMPERATURE DEPENDENCE OF THE 
RESISTANCE 



As temperature increases, the Fermi-Dirac distribution 
is broadened. Consequently the conductance of the of the 
transparent quantum point contacts (e c < fj) decreases 
exponentially (towards half its value), while that of the 
insulating ones increases. Thus we expect to see rather 
dramatic effects as a function of temperature. This is 
indeed depicted in Fig. 2. In (a) we plot the prediction 
of the model and in (b) the experimental data pj. As 
temperature is lowered, systems with slightly different re- 
sistance at high temperatures will diverge exponentially 
with decreasing temperatures. The resistance of systems 
on the metallic side (n > n c ) will saturate at zero tem- 
perature, while that of insulating samples will diverge, 
in agreement with the general shape of the experimental 
curves. Note that there is an upward turn even on the 
metallic side of the transition. In fact, close to the tran- 
sition, on the metallic side, as temperature decreases, the 
conductance of the insulating part decreases significantly, 
and its contribution to the total conductance is dramat- 
ically reduced. Since the critical percolation cluster is 
very ramified (in fact of fractal dimension), the contri- 
bution of the insulating part of the system dominates 
at high temperatures, and the increase of its resistance 
with decreasing temperature leads to an increase of the 
resistance as temperature is lowered, even on the metal- 
lic side. At low enough temperatures, however, when the 
resistance of the insulating part of the network becomes 
high enough, its contribution to the total conductance 
becomes negligible. Then the total resistance is domi- 
nated by the percolating conducting network, and thus 
the overall resistance will decrease with decreasing tem- 
perature. This leads to a nonmonotonic temperature de- 
pendence on the metallic side of the critical point, which 
can be clearly observed both in the model (Fig. 2c) and 
in the data (Fig. 2d). 

The fact that only deep in the metallic regime, the 
overall resistance increases with increasing temperature, 
suggests that the density at which the resistance is ap- 
proximately temperature independent is not the true crit- 
ical point, but rather deeper on the metallic side. This is 
clearly seen in Fig. 3b, where one can see a point where all 
the low-temperature curves nearly cross, well inside the 
metallic regime. The above discussion suggests that one 
should be cautious in associating the critical point with 
the experimentally observed "temperature-independent" 
point (Fig. 3a), as done routinely in the experiment in- 
terpretations. 
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Fig. 2. Temperature dependence of the resistance for sys- 
tems of different densities, as obtained by the model (a) and 
compared to the experimental data of Ref. |^| (b). The crit- 
ical line is denoted by the bold curve. All curves below the 
critical line saturate at zero temperature, while above it the 
resistance diverges. For systems close to the transition on 
the metallic side, the resistance is a nonmonotonic function 
of temperature, as seen both in the model (c) and in the data 
(d). 
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Fig. 3. Comparison of the density dependence of the con- 
ductance between the data of Q| (a) and the model (b), for 
several temperatures. The density at which the theoretical 
curves seem to cross each other, is well above the true critical 
point. 

Lastly, The resistance at the metallic regime is given by 
some geometrical factor times the inverse of a m (Eq. |^) , 
which naturally gives the observed exponential temper- 
ature dependence observed experimentally. The high- 
temperature resistance of the critical density network is 
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naturally around h/e 2 , the only resistance scale in this 
model. 



V. PARALLEL MAGNETIC FIELDS 



The effect of a parallel magnetic field on the overall 
conductance is determined by the way it affects the in- 
dividual points contacts. The effect of a parallel field 
on transport through a single QPC has been studied in 
detail 24 1. These experimental and theoretical studies 
demonstrated that the threshold density where the QPC 
opens up increases parabolically with the in-plane mag- 
netic field. This effect was attributed to coupling of the 
in-plane motion to the strong confinement in the vertical 
direction, leading to an increase in the confining enrgy. 
Writing, for simplicity, the three dimensional Hamilto- 
nian that describes free motion in two-dimensions and a 
harmonic confining potential in the third (z) direction, 
with a magnetic field pointing in the ir-direction 
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it is straightforward to se e that the bottom of the 2d band 
shifts from JmjJq/2 to H^Juj 2 + uj 2 /2, with lo c = eB/mc, 
leading to a corresponding decrease in the kinetic energy 
of all electrons. Thus the effective critical Fermi energy, 
or density, becomes larger. In other words, for a given 
density or chemical potential, if the system at zero field 
is on the metallic side, i.e. if the Fermi momentum is 
above the critical momentum (or kinetic energy) allow- 
ing percolation through the system, a parallel field will 
lower that energy towards the critical energy, eventually 
crossing the critical point and leading to an insulating 
behavior. Fig. 4 depicts the experimental data of p5[ , 
and the corresponding predictions of the model. As ex- 
pected, as magnetic field increases, the system gradually 
crosses over from a metallic to an insulating behavior. 

The above discussion allows a quantitative prediction 
of the effective critical energy in a parallel field , 



e c (H) = e c (H = 0) + h (^f^+^l - cj ) /2. 



(6) 



At zero temperature, when the resistance on the in- 
sulating side is infinite, we expect the resistance on the 
metallic side /i > e c (H = 0) to diverge with increasing 
field, 



R(H)~(u.-e c (H)Y 



(7) 
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Fig. 4. Comparison between the experimental data j25| and 
the predictions of the model, demonstrating that a parallel 
magnetic field causes the metallic system to cross over grad- 
ually to the insulating regime. 

For a finite temperature this divergence is cut off by 
the finite resistance of the insulator, and the magnetic 
field dependence changes as one crosses into the insulat- 
ing side. This behavior is clearly seen in the experimen- 
tal data of P5| ( very similar data was also reported by 
Mertes et al. p6[). Fig. 5 depicts the experimental data 
|p5j for the magnetic field dependence of the resistance, 
compared to what is expected from the model. The dif- 
ference in behavior between the metallic and the insulat- 
ing regimes is clear. On the metallic side we see that the 
resistance increases rapidly as the magnetic field brings 
the critical point closer to the chemical potential and 
then an abrupt change of behavior as the system enters 
the insulating regime. If the system was on the insulating 
side to begin with, then the magnetic field dependence 
of the resistance is similar to what is usually seen in sys- 
tems where transport is via variable-range hopping p7[ . 
In these system the positive magnetoresistance is due to 
spin polarization pq . The magnetic field dependence on 
the insulating side depicted in Fig. [v| is assumed here 
to be due to that process. In this regime, however, the 
magnetic field at which there is a marked change in be- 
havior is spin-related and should depend only weakly on 
density. A recent experimental investigation of the mag- 
netoresistance in the insulating side |29| indeed supports 
this mechanism. 
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Fig. 5. Comparison of the experimentally measured resis- 
tance, as a function of parallel magnetic field to the model 
predictions. On the metallic side, the magnetic field shift the 
critical point towards the chemical potential, leading to a di- 
vergence in the resistance, which is cut off by finite tempera- 
ture. 

As was mentioned above, the critical point in the den- 
sity - magnetic field plane shifts towards hi ghe r densities 
with increasing magnetic field. Yoon et al. [E5j have also 
measured the dependence of the critical field on density, 
which can be deduced from the inversion of Eq. (^J) , 

H c = m*c/e^&{H) 2 + 2hiu A(H), (8) 

where A(_ff) = e c (H) — e c (H = 0). The comparison of 
the prediction of this simple equation to the data is de- 
picted in Fig. 6 for two samples cut from the same wafer. 
The fitting parameters are the zero-field critical point, 



that can be read directly from the data, the gate capaci- 
tance - the rate at which the Fermi energy changes with 
density, and the confining energy in the perpendicular 
direction, hujQ. It is encouraging to note that while the 
critical energy, which is determined by the disorder real- 
ization, and the gate capacitance, which is determined by 
the geometry, arc different for the two samples, both sets 
of data can be fitted by the same value of the perpendic- 
ular confining energy, which ought to be the same for the 
two samples, and turns out to be hu>o ~ 0.8meV, leading 
to an extension of the wavefunction in the perpendicular 
direction of the order of 11 nm, similar to the value used 
by the authors of Ref. |u| to fit their experimental data. 
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Fig. 6. Comparison of the measured density dependence of 
the critical magnetic field [^5j (circles) to the prediction of the 
theory (Eq.^, solid lines). The two data sets are two different 
samples cut from the same wafer, and can be fitted using the 
same value of the perpendicular confining energy. 

Interestingly, it seems that the effects of parallel fields 
can be understood without employing the electron spin. 
As a parallel field will also reduce the conductance of 
some of the point contacts from 2e 2 /h to e 2 /h, the Zee- 
man effect will also increase the system's resistance. It 
should be noted that while the coupling of the in-plane 
motion to the confining potential in the perpendicular 
direction was also considered by Das Sarma and Hwang 
|p0|| , the magnetoresistance predicted here, in contrast 
with Ref. pc| , should not exhibit any anisotropy. The 
reason is that the direction of transport through the 
quantum point contacts is expected to be random, with 
no preferred direction. 

VI. PERPENDICULAR MAGNETIC FIELDS 

While the longitudinal resistance depends exponen- 
tially on temperature, the weak field Hall resistance is 
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practically independent of temperature |31|. Such an ob- 
servation might be hard to account for in theories that 
argue for new non-Fermi liquid-like behavior, but it is 
trivial in the present model - the critical exponent for the 
Hall coefficient in a two dimensional percolation problem 
is exactly zero p^] , and thus the Hall coefficient should 
display no critical behavior at the critical point. This 
prediction was indeed confirmed in classical percolation 
experiments pa], a nd is very similar to that observed by 
Pudalov et al. [|3lH . 
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Fig. 7. The critical density - the number of electrons in 
the puddle, so that the topmost energy will allow transport 
through the point contact - as a function of magnetic field, 
in the presence of a finite disorder. The continuous curve is 
an averaged fit through the (necessarily integer) data points. 
Top-right inset: the corresponding experimental data ifjjf . 
Bottom- left inset: results of effective medium theory 

The situation in larger perpendicular magnetic fields 
is more interesting, as quantum Hall (QH) states are 
formed. Transport through a single QPC in perpendic- 
ular field and the crossover between the zero field limit 
and the QH limit have been studied in detail J34| . As ex- 
pected, one finds that the critical energy oscillates with 
magnetic field due to the depopulation of Landau levels. 
In the present case, the oscillations are smoothed out by 
the disorder and by the averaging over many QPCs. Thus 
only the strongest oscillation, near v = 1, may survive, 
leading to a single dip in the critical density vs. magnetic 
field plot, as was observed experimentally. In order to al- 
low for the averaging procedure, one has to take the full 
conductance distribution into account, which is beyond 
the two-species scaling theory. For completeness I report 
here results of numerical calculations and effective 
medium theory [B5 36 1. In the numerical calculation I 



studied the energy levels of one puddle of electrons, which 
we modeled by a circular disk, in the presence of disorder 
p7|. In Fig. 7 I plot the "critical density" - the number 



of electrons that need to occupy the puddle, so that the 
energy of the highest-energy electron will be enough to 
transverse the QPC J3£|, equivalent in the bulk system 
to the critical density - as a function of magnetic field. 
Indeed a dip near v = 1 is clearly seen, with all other os- 
cillations smoothed out by the disorder. This curve has 
a strong resemblance to the experimental data |?9| (top- 
right inset). The results of the effective medium theory 
are depicted in the bottom- left inset, again demonstrat- 
ing a dip near v = 1. In addition, it is expected that as 
the magnetic field is lowered below the v — 1 minimum, 
more than one channel will transverse some QPCs, lead- 
ing to an increase in the critical conductance, as indeed 
reported experimentally. 



VII. REDUCTION FACTORS LARGER THAN 2 

In the degenerate electron limit (T <C Ep), the biggest 
reduction factor in the resistance on the metallic side, 
with decreasing temperature is a factor of 2. On the 
other hand, reduction factors close to an order of mag- 
nitude and even larger |4(| have been observed exper- 
imentally, especially in Silicon based samples . Here I 
show that when temperature becomes of the same order 
of magnitude of Ep, the reduction factor in the model 
can assume arbitrarily large factors ^] . 

The electron density is given by 



de 



1 + exp[(e - p)/T] 



Po Tlog [1 



■ exp(/VT)] 
(9) 



where po, assumed constant, is the electronic density of 
states (per energy and per volume), and /i, the chemical 
potential, is measured relative to the bottom of the band. 
Inverting the above equation, the chemical potential for 
a given density, is 

fi = Tlog [exp(n/p T) - 1] = Tlog [exp(E F /T) - 1] 

(10) 

where the Fermi energy is defined as the T — *• limit 
of the chemical potential. These textbook expressions 
demonstrate that while the Fermi energy varies linearly 
with the density, the chemical potential may be more 
sensitive to density variations in the nondegenerate limit 
T ~ Ep- Moreover, the chemical potential is now tem- 
perature dependent, and decreases with increasing tem- 
perature (see inset in Fig. 8). Substituting the above 
expression in the expression for the conductance through 
a single point contact (0) demonstrates that the conduc- 
tance can decrease arbitrarily with increasing tempera- 
ture (or, equivalently, that the resistance can decrease 
by an arbitrary factor with decreasing temperature). In 
Fig. 8 the temperature dependence of the conductance 
is plotted for two values of the Fermi energy, Ep. The 
curve for Ep T shows the expected behavior for the 
degenerate electron gas - the conductance decreases and 
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saturates at a value smaller from the zero temperature 
conductance by a factor of 2. On the other hand, in 
the nondegenerate regime, Ep ~ T, the conductance de- 
creases by a much larger factor. 
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Fig. 8. The temperature dependence of the conductance in 
the degenerate (kT <C Ep) and in the degenerate regime 
(kT ~ Ef). While in the nondegenerate regime, the con- 
ductance can decrease by a factor of two, it can decrease ar- 
bitrarily in the nondegenerate regime, due to the temperature 
dependence of the chemical potential (inset). 



VIII. CONCLUSIONS 

All the above results and discussion demonstrated that 
many of the experimental observations can be explained 
in the context of the simple semi-classical, noninteract- 
ing model introduced here. This is not to say that in- 
teractions and other effects are irrelevant. For exam- 
ple, the formation of the electron puddles may be dom- 
inated by interaction effects (see, e.g., Ref. fl3|| ), and 
dephasing is certainly dominated, at these low temper- 
atures, by electron-electron interactions. Other effects, 
including the energy dependence of the transmission co- 
efficient and the possibility of more than one channel 
through the QPCs, the role of interband-scattering ||] 
and temperature-dependent impurities ^ may also be 
important to understand quantitative aspects of the data. 
Nevertheless, the fact that several important aspects of 
the experimental data can be explained in the context of 
a simple model is quite encouraging. 

Some predictions made in where the model was 
presented, were already confirmed. The mechanism for 
the quenching of the metallic phase by a parallel mag- 
netic field was suggested there, and, as discussed above, 
agrees very well with recent experiments. Moreover, it 
was suggested that local measurements will be able to 
explore the percolative nature of the insulating phase. 



Indeed, Hani et al. |15| have used local probes to measure 
the change of the local chemical potential with density. 
While on the metallic side the signals from all probes were 
identical, and were accurately described by Hartree-Fock 
theory, these probes gave different signals on the insu- 
lating side, which the authors interpreted as a signature 
of a percolative phase. (An indirect experimental verifi- 
cation of the percolation process in the QH regime was 
already reported in [ff2|). As the metallic puddles can be 
thought of as quantum dots, one can use the abundant 
information about such structures p3| , to gain additional 
understanding of the characteristics of the puddles and 
the phase separation. Such local probes |Q can give a 
" smoking gun" verification of the picture presented here 
by looking for the periodic oscillations of the local chem- 
ical potential on the insulating side, due to depopulation 
of the Landau levels, as was observed in quantum dots 

To conclude, a semi-classical model, combining local 
quantum transport and global classical percolation was 
employed to explain the observed metal-insulator tran- 
sition. The model attributes the transition to the finite 
dephasing length in these temperatures. As temperature 
is lowered even further, and the dephasing length be- 
comes larger than the puddle size, quantum localization 
effects should kick in. Such weak-localization corrections 
in the metallic side were indeed observed experimentally, 
confirming the expectation that, if the dephasing length 
indeed diverges at zero temperature, these systems will 
eventually becomes Anderson insulators. 
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